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ABSTRACT 

We  discuss  the  well-posedness  of  the  model  problem 

ut  =  *(ux}x  °n  10,15  X  t0,T5#  T  >  0  '  (P> 

subject  to  given  Neumann  or  Dirichlet  boundary  conditions  at  x  =  0  and 

x  ■  1 ,  and  to  the  initial  condition  u( x,0 )  -  f (x) i  the  given  functions 
f  :  [0,1]  ♦  R,  #  :  R  ♦  R  are  assumed  to  be  smooth/  $(0)  *  0,  $  satisfies 
the  coercivity  assumption  K)  >  c?2,  for  some  constant  c  >  0  and  for 
€  e  R,  and  #  is  assumed  to  be  decreasing  on  an  interval  (a,b)  with 
a  >  0.  We  present  a  recent  nonuniqueness  result  in  the  special  case  when  # 
is  piecewise  linear  and  study  a  related  convexified  problem. 
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SIGNIFICANCE  AND  EXPLANATION 


The  equation  (P)  in  the  abstract  can  be  viewed  as  a  simple  model  for  non¬ 
linear  diffusion.  The  existing  mathematical  theory  requires  >  0.  However, 
in  one  space  dimension  the  laws  of  thermodynamics  merely  imply  that  the  graph 
of  $  lies  in  the  first  and  third  quadrant,  without  necessarily  requiring  4 
to  be  monotone  nondecreasing.  This  raises  the  natural  question  whether  the 

assumption  >  0  can  be  replaced  by  the  much  weaker  coercivity  condition 
2 

£♦(5)  >  c£  ,  c  >  0.  For  a  nonmonotone,  piecewise  linear,  coercive  #  it  was 
shown  in  MRC  TSR  #2354  (see  [5] )  that  the  initial  value  problem  for  (P)  has 
infinitely  many  solutions,  whenever  the  initial  data  f ' ( •)  reach  the  criti¬ 
cal  range  where  #* (•)  <  0.  These  solutions  u  of  (P)  have  the  property  that 
^  oscillates  between  regions  in  which  ♦' ( •)  >  0.  Although  CP)  is  evidently 
not  well-posed,  it  is  hoped  that  imposing  additional  physically  motivated 
assumptions  will  lead  to  a  natural  selection  and  a  well-posed  problem. 

In  this  report  we  first  review  known  (previously  unpublished)  a  priori 

estimates  for  (P),  and  we  then  give  a  simpler  construction  for  the  existence 

p*  / 


of  infinitely  many  solutions  of  (P)  for  a  piecewise  linear  £  as  suggested  by 
G.  Strang  We  then  investigate  further  the  qualitative  behavior  of 

solutions  of  (P).  Motivated  by  known  results  in  one  space  dimension  for  the 
steady  state,  non-elliptic  problem^  we  study  the  analogous 

convexified  problem  associated  with  (P).  Analytical  and  numerical  considera¬ 
tions  suggest  that  the  unique  solution  of  the  convexified  problem  (which  has  a 
monotone,  nondecreasing  constitutive  function)  can  be  interpreted  as  an 
average  of  solutions  of  CP) ,  whenever  the  data  f ' (  •)  reach  the  critical 
range. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  authors  of  this  report. 
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1 .  INTRODUCTION 

We  discuss  the  Initial  boundary  value  problem 

ut  *  on  l0*1!  *  [0»T] 

ux(0,t)  »  u^l.t)  -  0,  u(x,0 )  »  f(x)  , 

where  subscripts  denote  partial  derivatives,  under  the  principal  assumption: 
♦  i  R  ♦  R  smooth,  4(0)  »  0,  and  there  exists  a 
constant  c  >  0  such  that  54(C)  >  cC2 ,  (  e  1, 

and  there  exists  an  interval  (a,b)  with  a  >  0  (A) 

such  that  ♦’ (5)  <  0,  ?e  (a,b). 

A  model  case  is  sketched  in  Figure  1 . 


Figure  1 


Concerning  f  we  assume  throughout  that  f  :  [0,1]  ♦  R  is  as  smooth  as 
needed  and  satisfies  the  boundary  conditions.  Dirichlet  and  inhomogeneous 
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boundary  conditions  can  also  be  discussed,  including  the  theorems  in  Sections 
3  and  4. 

One  motivation  for  studying  (F)  is  the  theory  of  nonlinear  diffusion. 

Specifically  the  Clausius-Duhem  inequality  [3,  p.  79]  in  one  space  dimension 

can  be  shown  to  imply  that  the  flux  4(ux)  need  not  be  a  monotone  function  of 

the  gradient  of  temperature  ux;  this  is  consistent  with  our  assumption  (A) . 

If  $  is  strictly  monotone  increasing  (as  is  usually  assumed)  standard  theory 

guarantees  that  (P)  has  a  unique  solution  which  is,  roughly  speaking,  as 

smooth  as  the  function  In  particular  ( see  Figure  1 ) ,  this  is  true  in  the 

model  case  if  the  data  f  satisfy  f'(x)  <  a.  If,  however,  (f'(xQ))  <  0 

for  some  xQ  e  (0,1),  then  in  a  neighborhood  of  (xq,0),  (P)  behaves  like  a 

backward  parabolic  equation  which  is  not  well  posed.  In  particular,  if  in 
13  3  J 

Figure  1  ♦(5)”y£  -  —  t  +  2C»  then  (P)  cannot  have  a  solution  u  such 
that  ux  is  piecewise  continuous  on  [0,1]  x  [0,T]  for  any  T  >  0,  unless 
f  is  analytic.  The  basic  problem  to  be  discussed  here  is  whether  (P)  is 
well-posed  in  some  precise  sense,  whenever  the  data  f'(x)  fall  in  the 
critical  range  ((a,b)  in  the  model  case). 

The  a  priori  estimates  given  in  Section  2  suggest  that  (P)  may  have 

m 

solutions  u  with  e  L  ([0,1]  x  [0,T] )  for  some  T  >  0.  Such  solutions 
are  constructed  in  Section  3  for  piecewise  linear  $  satisfying  assumption 
(A) ;  however,  the  solutions  are  not  unique.  Motivated  by  known  results  for 
the  steady-state,  non-elliptic  problem  ( Wux)x  *0),  we  are  lead  to 
discussing  in  Section  4  the  "convexif ied"  problem  associated  with  (P),  with 
the  objective  (not  yet  established)  of  distinguishing  a  solution  of  (P)  with 
special  properties.  Other  open  questions  include  the  existence  of  solutions 
of  (p)  when  4  is  not  piecewise  linear,  and  a  justification  of  phase  changes 
suggested  by  numerical  experiments. 


2.  REMARKS  ON  A  PRIORI  ESTIMATES 

Consider  (P)  with  either  homogeneous  Neumann  or  Dirichlet  boundary 
conditions  at  x  *  0,1,  and  with  f  smooth.  Let  u  e  W1 ' 2  be  a  solution  of 
•  (P)  on  [0,1]  x  [o ,T]  for  some  T  >  0.  Standard  energy  methods  yield  the 

following  a  priori  information. 

a.  Multiplying  by  u  and  integrating  the  p.d.e.  over  [0,1]  *  [0,T]  gives 

1  1  1 

(i)  ess  sup  /  u2(x,t)dx  <  —  /  f2(x)dx,  and 

te[0,T]  0  0 


-3 


□  □ 


It  should  be  observed  that  If  +  Is  strictly  monotone  these  estimates 
imply  that  u  e  H  ,  and  standard  techniques  yield  existence,  without  applying 
e.g.  the  theory  of  maximum  monotone  operators  (in  particular,  this  is  true  if 
f'(x)  <  a,  see  Figure  1).  However,  one  cannot  do  this  if  t  is  not  monotone 
everywhere  (+'(f'(x))  <  0  for  some  x  e  R). 

A  more  subtle  estimate,  proved  by  J.  Bona,  L.  Wahlbin  and  the  second 
author,  is 

For  either  the  Dirichlet  or  Neumann  problem 

*U  1  <  1  sup  |+(f'(x))|  t 

xe[0,1] 

the  optimal  choice  of  the  coercivity  constant  c  ■»  inf  >  0. 

yen  y 

2k 

Proof:  Multiply  the  p.d.e.  by  +  ^ux^ut'  wher®  in  th®  exponent  k  is  a 
positive  integer,  and  integrate  with  respect  to  x  and  t.  Defining 


V(y)  *  /  ♦2k+1  ( £)d£  (>  0  by  hypothesis  on  +  for  y  C  R) 


one  obtains,  after  using  integration  by  parts  and  the  boundary  conditions: 


t  1 


f  f  *2k<\>\****  +  zrh  1  VVdx "  irh  [  Vf,(x))dx  *  (1) 


0  0 


From  the  coercivity  assumption 


2k+1 

fk(y>  >  2k  ~ 2  y2k+2  (y  e  and 


1  2k+1  1 

/  T  (u  ) dx  >  - --■  /  u?k+2dx 

i  kx  2k  +  2 


-4- 


Next  let  It-  »«•«-»-  sup  U(£'(x)|»  then 

xeto,i] 


/  »  ( f'(x))dx  < 


Dropping  the  positive  first  term  in  (1),  and  then  raising  both  sides  to  the 


(2k  +  2)"1  yields 


L  _1 _  2k+1 

<  »2~2  »2k‘2. » <  *  ‘ T  • 


(2k  +  2) 

finally «  letting  k  ♦  -  we  obtain 

lul  <  -  ,  0  <  t  <  T  • 

X  -  c 

d.  Mgisgainimia  principle.  Consider  (P)  in  the  model  case  of  Figure  1 

with  homogeneous  Neumann  boundary  conditions*  Let  u  fi  W  be  a  solution 

on  [0,1]  x  [0,T]  for  some  T  >  0.  Then  for  0  <  t  <  T: 

mints,  min  f'(x))  <  u*(x,t)  <  maxtb,  max  f'(x))  . 

0<x<1  0<x<1 

For  classical  solutions  the  result  can  be  proved  by  a  standard  comparison 
method.  We  sketch  a  more  general  approach.  Let  v  *  ux.  Then  v  is  a 

solution  of  the  Dirichlet  problem 

v  *  *(v)  on  [0,1]  x  [0»T] 

t  "  XX  (2) 

v(x,0)  -  f'(x)i  v(0,t)  -  v(1,t)  -  0  . 

Let  0  <  a'  <  a,  b'  >  b  >  0  (see  Figure  1).  Define  g  :  *  ♦  *  such  that 

s 

g. (.)  .  o  on  (a* ,b* )  and  g'(*)  >  0  otherwise.  Let  G(s)  -  /  g(  C)dC. 
Multiply  (1)  by  g(v)  and  integrate  over  [0,1]  x  [0,t],  0  <  t  <  T. 
Integration  by  parts  and  the  boundary  conditions  yield 


1  t  1  2  i 

/  G(v(x,t)  )dx  +  /  /  (v)g’  ( v)v*dxdt  -  /  G(fMx))dx  . 


-5- 


By  definition  of  g  and  assumptions  on  4  the  double  integral  is  positive , 


so  that 


1  1 

/  G(v(x,t))dx  <  /  G(f'(x))dx,  0  <  t  <  T  . 

0  0 


The  successive  specific  choices 

g(.)  .  (.  -K)  ,  k  -  max(b' ,  max  f'(x))  , 
r  0<x<1 


and 


g(  •)  *  - (jc-  •  ).,  k  ■  min(a' ,  min  f'  (x)) 

0<x«1 


yield  G(v(x,0 ) )  -  G(f'(x))  -  0.  Therefore,  respectively. 


/  (v<x,t)  -  k)2dx  <0,  0  <  t  <  T  , 


and 


\  /  <k  -  v(x,t))2dx  <0,  0  <  t  <  T  . 

0 


Thus 

k  <  v(x,.t)  <  k,  0  <  t  <  T  . 

Letting  a'  ♦  a  and  b'  ♦  b  yields  the  result.  It  should  be  noted  that  this 
stronger  result,  proved  under  the  stronger  assumption  u  e  W2'2  (rather 
than  w1'2),  is  not  implied  by  the  L**  estimate  in  2c. 


3.  EXISTENCE  AND  NONUNIQUENESS  IN  A  SPECIAL  CASE 

For  a  piecewise  linear  constitutive  function  $  it  was  shown  in  [5]  that 
the  problem  (p)  has  infinitely  many  solutions  if 

{x  -  ♦'(f’(x))  <  0)  *  4  :  (3) 

Theorem.  Let  $  be  piecewise  linear  with  ((-•,a)  U  ( b , ») )  >  0  and 
t'((a,b))  <  0  (see  Figure  1)  and  assume  that  f'(x)  satisfies  (3).  There 
exists  T  >  0  such  that  (P)  has  infinitely  many  solutions  u  with 
HSlder  continuous ,  and  ux6L,utei^(p<  •)  on  the  domain 
[0,1]  x  [0,T]  •  Moreover,  u^Xjt)  0  [a,b]  for  (x,t)  e  [0,1]  x  (0,T]. 

We  outline  the  proof  of  the  Theorem  in  the  case  when  the  increasing  parts 
of  t  Are  parallel,  i.e.  $  is  of  the  form: 


Figure  2 

with  +'<K\[a,b])  -  1.  For  simplicity  we  also  assume  that  f'(x)  <  2,  x  e  [0,1]. 
The  key  to  the  proof  is  the  construction  of  a  function  w  which  represents 
the  oscillating  part  of  a  solution  u  and  reduces  the  problem  ( 1 )  to  an 
inhomogeneous  heat  equation  for  the  smooth  part  v  of  the  solution  u  -  v  +  w. 

We  describe  an  idea  of  6.  Strang  [6]  for  constructing  w  which 
simplifies  the  original  argument  in  [5].  Consider  the  following  (infinite) 
triangulation  of  [0,1] 2 
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which  is  determined  by  the  points 

zjk  (k2  ^'2  j)'  k  "  j  e  *  . 

Set  h(0)  =  0,  h* (x)  -  max(0,f’(x)  -  1).  Define  w  as  the  piecewise  linear 
function  with  respect  to  this  partition  which  interpolates  the  data 
Wjk  s-  max{v2  ^  s  V2  ^  <  h(k2  "*),  v  e  *} 
at  the  points  z^. 

By  the  definition  of  w^, 

Wjk  <  h(k2~^)  <  w.^  +  2~^  , 
and  since  0  <  h’  <  1  we  have 

wj(k+1)  "  wjk  e  *0'2  ^  • 

One  then  easily  verifies  that 

(i)  wx  s  [0,1]  x  [0 ,T]  {0,1} 

(ii)  wt  e  l" 

(iii)  lim  lw( «,t)  -  hl^  -  0(4) 
t+0 

(iv)  f(x)  <  1  -  2"3  — >  wx(x,t)  -  0  for  t  <  2~*  . 

Property  (iv)  means  that  for  small  t  wx( »,t)  is  supported  in  a  neighborhood 
of  the  set  {x  s  f(x)  >  1}.  Let  v  be  the  solution  of  the  problem 


-8- 


(5) 


vx(0,t)  »  vx(1,t)  »  0 
•  v(x,0)  -  f(x)  -  h(x)  . 

We  claim  that  u  ■  v  +  w  is  a  solution  of  problem  (P). 

Since  wfc  6  lT ,  vx  is  Hfilder  continuous  and  vfc  e  iP  for  any  p  <  “  by 
standard  estimates  for  the  heat  equation.  Put  ic  min(a  -  1,2  -  b).  By  the 
continuity  of  vx  and  (4{iv))  there  exists  T  >  0  such  that 
vx(x,t)  <  1  +  x,  (x,t)  e  [0,1]  x  [0,T]  , 

and 

(supp  wx)  n  [0,1]  x  [0 ,T]  C 

{( x, t)  :  t  e  [0,T] ,vx(x,t)  >  1  -  x)  . 

In  view  of  the  relation 

+(£  +  1)  -  *(C>,  5  e  [1  -  x,1  +  «3  ,  (6) 

l 

it  follows  that 

v.  +  w.  «  v,  **  $(v  )  “  $(v  +  w  )  . 

t  t  XX  T  X  X  X  XX 

Clearly,  the  function  w  is  not  uniquely  determined.  For  any  X  >  0, 
w^(x,t)  »•  v(x,Xt)  satisfies  (4).  Since  vx  is  continuous  the  solutions 
”  vx  +  ar®  distinguished  by  their  different  discontinuity  patterns. 
Thus  there  is  a  continuum  of  solutions  to  the  problem  (P). 

If  the  monotone  increasing  parts  of  $  have  different  slopes  the  proof 
of  the  Theorem  is  considerably  more  complicated  (see  [5]).  The  reason  is  that 
instead  of  (6)  we  have  the  relation 

♦U  +  (A5  +  B))  -  ♦(£)  ,  (7) 

where  the  shift  is  no  longer  independent  of  €•  Nevertheless,  the  proof  is 
based  on  the  construction  of  an  auxiliary  function  w  which,  however,  depends 
on  v  and  the  equation  for  the  smooth  part  v  becomes  quasilinear.  For  a 
general  smooth  nonmonotone  $  the  existence  question  is  not  settled. 
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fa 


In  spite  of  nonuniqueness,  numerical  approximations  to  the  problem  (P) 
are  stable;  the  x-derivatives  of  the  numerical  solution  have  the  oscillatory 
behavior  asserted  by  the  Theorem. 

We  used  a  Crank-Nicholson  scheme 

uj  uj  “  to  (Duj-1  +  "“j-l*  )  (8) 

where  u”  approximates  u( j Ax, nAt)  and  DUj  “  u^+1  -  Uj  denotes  the  forward 

difference.  For  $  in  Figure  2  with  a  =  1.25,  b  -  1.75  and  the  model 

initial  data  f(x)  =  3x2  -  2x3  the  numerical  solutions  show  the  following 

typical  behavior.  For  small  t,  uAt  <  0.1,  Du^  oscillates  over  the  interval 

where  (f ' )  is  negative  as  shown  in  figure  (4)  below 


—5  2  -.7 

where  Ax  »  2  ,  At  -  (Ax)  ,  nAt  “  2  .  The  oscillations  are  very  regular. 

Du"  alternates  between  the  values  1  and  2  which,  in  agreement  with  the 
Theorem,  differ  by  1  and  do  not  lie  in  the  interval  [1.25,  1.75)  where  ♦' 
is  negative.  Moreover,  $(Du”)  is  (numerically)  continuous. 

As  t  ■  nAt  increases  the  oscillations  gradually  disappear  and  u" 
approaches  the  limiting  state  ux  *  0.  Similar  numerical  approximations  have 
been  independently  studied  by  G.  Strang  and  A.  Naby  [6].  It  has  not  been 
established  that  the  numerical  approximations  converge  and  thereby  possibly 
single  out  a  distinguished  solution  of  the  problem  (p). 
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4.  A  RELATED  CONVEX! PIED  PROBLEM 


It  is  known  [1]  that  in  one  space  dimension  the  solution  of  the  problem 


0  *  *(ux)x'  x  e  10,11 
u(0)  -  0,  u(1)  *  d,  d  a  constant  , 


(9) 


can  be  obtained  by  minimizing  the  functional 


1 

J(u)  «  /  4(u  (x))dx,  4' {  •)  ”  4<  •)  , 

0 


over  all  functions  u  6  W1 ,2(  [0,1] )  satisfying  the  given  boundary  conditions, 
under  suitable  assumptions  on  4  (these  do  not  require  4  convex).  If  4 
satisfies  assumptions  (A),  (9)  is  the  steady  state  problem  associated  with 
(P).  By  a  result  of  Ekeland  and  Tfaam  [4]  such  a  solution  u  of  the 
variational  problem  also  minimises  the  functional 


1 

5(u)  -  /  4(u  (x))dx  , 
0  x 


where  4  is  the  convexification  of  4,  and  u  solves  the  convexified 
problem: 


0  -  4<«x>x*  *  «  I®'1} 

.9 


(10) 


u(0)  -  0,  u(1)  -  d,  4(  •)  *  4  (  •)  f 

the  function  +(•)  is  necessarily  nondecreasing  on  R  (if  4  has  the  graph 
in  Figure  2,  Figure  5  below  is  the  graph  of  4).  An  analogous  variational  prin¬ 
ciple  has  been  established  for  problem  (p),  but  only  for  4  nonotone  (see  [2]). 
Motivated  by  these  results  we  consider  the  convexified  problem  associated 


with  (P): 

Ut  -  4(ux>x»  e  CO , t ]  x  [0,T]  , 

u^O.t)  «  ux(1,t)  -  0,  u(x,0)  -  f(x)  . 


(5) 
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By  the  Theorem  below,  solutions  of  (P)  and  (?)  show  a  qualitatively  different 

/ 

behavior  and,  at  first  sight,  there  is  ho  obvious  connection  between  them  when 

/ 

the  data  f'  is  smooth  and  falls  in  tile  critical  range  (a,b).  Nevertheless, 
it  appears  plausible  that  as  t  ♦  »  problems  (P)  and  (?)  approach  the  same 
steady  state. 

For  simplicity  we  discuss  (?)  only  for  the  piecewise  linear  constitutive 
function  considered  in  Section  3.  For  $  in  Figure  2,  is  of  the  form 


Figure  S 


The  graph  of  4  is  familiar  from  the  Stefan  problem.  In  fact,  if  the  initial 
data  satisfy 

f'(x)  t  (1,2),  x  e  [0,1)  , 
t(f')  continuous  , 

(?)  describes  a  Stefan  problem.  If,  however,  f  is  smooth  which  we  assume 
throughout  this  section,  the  evolution  of  the  free  boundaries  resulting  from 
(?)  is  different.  For  simplicity  we  assume  that 

f*  (x)  <2,  x  e  [0,1)  , 

{x  i  f' (x)  >  1}  -  (rQ,Sg),  0  <  rQ  <  s0  <  1  ,  (11) 

f"(r„)  >  0,  f"(s0)  <  0  . 

Theorem.  Problem  (?)  with  $  given  in  Figure  5  has  a  unique  solution  u 
on  [0,1]  x  R+  with 

lu  I  <  If' I 

X  »  • 


-12- 


Iu„l  <  !♦(«•) '  »  • 

If  f  satisfies  assumptions  (11)#  the  region  fl  s “  {(x,t)  :  ux(x,t)  >1} 
where  “0  is  bounded  by  two  curves  r,s  6  ^  ( Hdlder  class)  connecting 

the  points  (r0,0) #<s0,0)  with  the  point  (x*,T*>  (Figure  6).  Moreover, 

T*  <  If'  .  (12) 


Note  that  u  satisfies  the  heat  equation  on  {0,1)  *1^X0  and  the  regularity 
assertions  in  particular  describe  the  behavior  of  u  in  a  neighborhood  of  the 
points  (r0,0),(s0,0). 

The  regularity  for  u,r,s  is  proved  by  using  the  implicit  semidiscrete 
approximation 

u(x,(n+1 ) At)  -  u(x,nAt)  •  Atf(u  (x,(n+1 ) At) )  ,  x  e  (0,1)  , 
of  the  equation  ut  -  ♦*ux)x*  Prove  the  inequality  (12).  Since  u  is 
continuous  and  u(x,t)  ”  u(x,0)  on  (2  we  have 

u(r(t),t)  -  f (r(t) )  .  (13) 

Also,  note  that 

ux(r(t)",t)  -  1  .  (14) 
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Integrating  the  equation  u.  -  $(u  )  over  the  domain 

^  XX 

Qy  «»  {(x,t)  t  0  <  x  <  r(t),  t  e  (0,T#)}  we  obtain 
0  ■  //  (u^  -  u^ldxdt  - 


-  /  f  (x)dx  -  / 

0  0 


r(T#) 

f(r(t) )r(t)dt  +  / 

0 


u(x,T*)dx 


/  1  dt  . 

0 


Using  (13)  and  lu^ <  If* 1^  for  the  estimate,  it  follows  that 

r(T*)  1  r(T*> 

T  *  **  J  (u(x,T#)  -  f(x))dx  <  /  /  2  If  *  l^dydx  <  If*  . 

0  Ox® 


Numerical  computations  show  (Figure  6  is  a  particular  example  with 
2  3 

f(x)  «  3x  -  2x)  that  the  curves  r,s  touch  the  x-axis.  Formally,  this  can 
be  explained  as  follows.  Differentiating  equation  (13)  with  respect  to  t 
»nd  using  (14)  we  have 


Since  f'(r(0))  - 


where  g(£  -  r0)  : 
continuous  on  Q 


r(t) (f ' (r(t) )  -  1)  *  ut(r(t),t)  . 
f'(r0)  «  1  we  can  rewrite  this  in  the  form 
r(t) (r(t)  -  rQ)  ■  u^rtt)  ,t)/g(r(t)  -  r0) 

«*<5)  -  f’(r0) 

“  - : -  .  If  we  assume  that  u.  ,u>v 

*  "  ro 

it  follows  that 


are 


r 

r(t)  -  rQ  +  /5t  +  o(^t)  (t  ♦  0)  ,  (15) 

which  is  consistent  with  regularity  of  r  established  in  the  theorem  above. 

We  have  not  yet  succeeded  in  proving  (15).  The  difficulty  is  that  the 
free  boundary  x  ■  r(t)  is  characteristic  for  (?)  at  the  point  (rQ,0). 
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